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Abstract 



Ventricular fibrillation (VF), the major reason behind sudden cardiac death, 

is turbulent cardiac electrical activity in which rapid, irregular disturbances 

in the spatiotemporal electrical activation of the heart makes it incapable of 

any concerted pumping action. Methods of controlling VF include electrical 

defibrillation as well as injected medication. Both methods yield results that 

are subject to chance. Electrical defibrillation, though widely used, involves 

subjecting the whole heart to massive, and often counterproductive, electrical 

shocks. We propose a defibrillation method that uses a very-low-amplitude 

shock (of order mV) applied for a brief duration (of order 100 ms) and over 

a coarse mesh of lines on our model ventricle. Our proposal is based on a 

detailed numerical study of a model for ventricular fibrillation ||T|. 
PACS numbers: 87.19.Hh, 05.45. Gg, 05.45. Jn, 82.40.Ck 
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Typeset using REVT^ 



Ventricular fibrillation (VF), the leading cause of sudden cardiac death, is responsible for 
about one out of every six deaths in the USA In the absence of any attempt at 

defibrillation, VF leads to death in a few minutes. The importance of this problem can 
hardly be overemphasized, so it is not surprising that many studies of VF have been carried 
out in various mammalian hearts and mathematical models 0]. Recent experiments 
have shown that, though VF is turbulent cardiac activity, it is associated with the 
formation, and subsequent breakup, of electrophysiological structures that emit spiral waves; 
these are referred to as rotors or simply spirals. Such structures and their breakup have 
also been obtained in a set of partial differential equations of the Fitzhugh-Nagumo type; 
this set has been proposed as a model for VF |l|]. We show that the spiral turbulence 
associated with VF in this model arises because of spatiotemporal chaos; this is similar to 
the spatiotemporal chaos in a related model for the catalysis of CO on Pt(llO) which we 
have studied elsewhere P,^. Our defibrillation scheme for the VF model of Ref. [|l| relies on 
the understanding we have developed for spatiotemporally chaotic states in the CO catalysis 
model. 

Electrical defibrillation entails the application of electrical jolts to the fibrillating heart to 
make it start beating normally again. Such defibrillation works only about two-thirds of the 
time and often damages the heart in the process of reviving it . In external defibrillation, 
electrical shocks 5 kV) are applied across the patient's chest; they depolarise all heart cells 
simultaneously and essentially reset the pacemaking nodes of the heart |TT|]. Slightly lower 
voltages (~ 600 V) suffice in open-heart conditions, i.e., when the shock is applied directly 
on the heart's surface [l^. In our study of the VF model IQ we achieve defibrillation by 
using very-low- amplitude pulses (of order mV) applied for a 6ne/ duration (of order 100 ms) 
and over a coarse mesh of lines on our model ventricle. Thus, if our defibrillation scheme 
can be realised in an internal defibrillator, it would constitute a significant advance in these 
devices. 

Before proceeding to our quantitative results it is useful to define the model for ventricular 
fibrillation [||]. For simplicity we use the model for isotropic cardiac tissue; in this case the 
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equations governing the excitability e and recovery g variables are 

de/dt = \/'e-fie)-g, 

dg/dt = e{e,g){ke- g). 
The function /(e), which specifies fast processes (e.g., the initiation of the action potential) 
is piecewise linear: /(e) = Cie, for e < ei, /(e) = —C2e + a, for ei < e < 62, and 
/(e) = C^i^e — 1), for e > 62- We use the physically appropriate parameter values given in 
Ref. Il, namely, ei = 0.0026, 62 = 0.837, Ci = 20, C2 = 3, C3 = 15, a = 0.06 and k = 3. 
The function e{e,g), which determines the dynamics of the recovery variable, is €{e,g) = ei 
for e < 62, e(e, (7) = €2 for e > 62, and e{e,g) = for e < ei and g < gi with gi = 1.8, 
ei = 1/75, 62 = 1.0, and £3 = 0.3 (which lies in the range 2 > 63 > 0.1 suggested in Ref. |jl|). 
We solve model (1) by using a forward-Euler integration scheme. We discretise our system 
on a grid of points in space with spacing 6x = 0.5 dimensionless units and use the standard 



five- and seven-point difference stencils |T3| for the Laplacians in spatial dimensions d = 2 
and 3, respectively. Our spatial grid consists of a square lattice with L x L points or cubic 
lattices with L x L x Lz points; in our studies we have used values of L ranging from 128 
to 512 and 2 < Lz < 16. Our time step is 6t = 0.022 dimensionless units. As in Ref. [|l|, 
we define dimensioned time T to be 5 ms times dimensionless time and 1 spatial unit to be 

1 mm, such that the period and wavelength of a spiral wave are approximately 120ms and 
32.5 mm, respectively. The dimensioned value of the conductivity constant in model (1) is 

2 cm^/s Our initial condition is a broken plane wave which we allow to evolve into a 
state displaying spiral turbulence; we control this eventually by our defibrillation scheme. 
On the edges of our simulation region we use no-fiux (Neumann) boundary conditions since 
the ventricles are electrically insulated from the atria. We impose these boundary conditions 
numerically by adding an extra layer of grid points on each side of our simulation grid and 
requiring the values of e be equal pointwise to their values on the points one layer within 
the boundary. 

Before constructing an efficient defibrillation scheme for VF in model (1) it is important 
to appreciate the following points: (a) Ventricular fibrillation in this model arises because 



the system evolves to a state in which large spirals break down [|l],H. (b) This state is a 
long-lived transient whose lifetime tl increases rapidly with L, the linear size of our system 
(e.g., for d = 2, ~ 850 ms for L = 100 whereas tl — 3200 ms for L = 128). This 
property of model (1) is in qualitative accord with the experimental finding that the hearts 
of small mammals are less prone to fibrillation than those of large mammals p!^JT5|]. For 



time t ^ Tl, a quiescent state with e = g = is obtained, (c) In systems with L > 128, 
Tl is sufficiently long that a nonequilibrium statistical steady state is established. This 
state displays spatiotemporal chaos. For example, we find that there are several positive 
Lyapunov exponents Aj (averages for Aj are performed for tq < t < t^, where tq is the time 
of decay for initial transients |l^); the number of positive Aj increases with L (e.g., the 
Kaplan- Yorke dimension Dky increases from 7 to ~ 35 as L increases from 128 to 256). 



(d) Model (1) is akin to a model for the catalysis of CO on Pt(llO) in so far as both show 
spiral breakup and spatiotemporal chaos [l|,§] . Our recent studies of the CO catalysis model 
have shown that Neumann or no-fiux boundary conditions tend to absorb spiral defects 
and, indeed, the spirals do not last for appreciable periods of time on small systems. 
Given that the long transient which leads to VF in model (1) is spatiotemporally chaotic, we 
might guess that the fields e and g have to be controlled globally to achieve defibrillation. In 



fact, some earlier studies |]T8[ of spiral breakup in models for ventricular fibrillation have used 
global control. Here we show that a judicious choice of control points (on a mesh specified 
below) leads to an efficient defibrillation scheme for model (1). For d = 2, we divide our 
simulation domain (of size Lx L) into smaller blocks and choose the mesh size such that 
it effectively suppresses the formation of spirals. For d = 3, we use the same control mesh 
but only on one of the square faces of the L x L x Lz simulation box. In our defibrillation 
scheme we apply a pulse to the e field on a mesh composed of lines of width 36x. A network 
of such lines is used to divide the region of simulation into square blocks whose length in 
each direction is fixed at a constant value L/K for the duration of control. (The blocks 
adjacent to the boundaries can turn be out to be rectangular). The essential point here is 
that, if a pulse is apphed to the e field at all points along the mesh boundaries for a time Tc, 



then it effectively simulates Neumann boundary conditions (for the block bounded by the 
mesh) in so far as it absorbs spirals formed inside this block (just like Neumann boundary 
conditions absorb spiral defects in the CO catalysis model |^). Note that is not large at 
all since the individual blocks into which the mesh divides our system are of a linear size 
L/K which is so small that it does not sustain long, spatiotemporally chaotic transients. 
Nor does which is related to the mesh density, have to be very large since the transient 
lifetime, r^,, decreases rapidly with decreasing L. We find that, for (i = 2, L = 128, K = 2 
and Tc = 41.2 ms is required for defibrillation. In Fig. 2 we illustrate such defibrillation 
with Tc = 44 ms. For d = 2, L = 512, K = 8 and = 704 ms suffices. Finally we show that 
a slight modification of our defibrillation scheme also works for d = 3 (Fig. 3). 
The efficiency of our defibrillation scheme is fairly insensitive to the height of the pulse we 
apply to the e field along our control mesh so long as this height is above a threshold. To 
obtain the value of e in mV units we have scaled the peak amplitude of a spike in the e 
field (which has an amplitude of 0.9 in dimensionless units) to be equal to 110 mV. The 
latter is a representative average value of the peak voltage of an electrical wave in the heart 



1^ . With this voltage scaling, dimensioned excitability is computed as 110/0.9 ~ 122.22 
mV times dimensionless excitability. We find, e.g., that, for L = 128, the smallest pulse 
which yields defibrillation is 57.3 mV/ms for the parameter values we use; however, we have 
checked that even stronger pulses (e.g., 278.3 mV/ms) also lead to defibrillation. We use a 
capacitance density of 1 yuF/cm^ |]2T|, which then yields a current density of 57.3 yuA/cm^. 
Note that this threshold is of the order of the smallest potential (~ 22.18 /iA/cm^) required 
to trigger an action-potential spike in model (1) with g = 0; the nuUclines for model (1), in 
the absence of the Laplacian, are such that this smallest potential increases with increasing 
g. This is physically consistent with the increase in the refractory nature of the heart tissue 
with increasing g. 

We have checked that (a) small, local deformations of our control mesh or (b) the angle of 
the mesh axes with the boundaries do not affect the efficiency of our defibrillation scheme. 
Furthermore the application of our control pulse on the control mesh does not lead to an 



instability of the quiescent state; thus it cannot inadvertently promote VF. We have checked 
specifically that, with e = 0, = as the initial condition at all spatial points, a wave of 
activation travels across the system when the control pulse is initiated; this travels quickly 
(~ 200 ms for L = 128) to the boundary where it is absorbed and quiescence is restored. Our 
defibrillation method also works if g is stimulated instead of e. This can be implemented by 
pharmaceutical means in an actual heart. 

Our two-dimensional defibrillation scheme above applies without any change to thin slices of 
cardiac tissue. However, it is important to investigate whether it can be extended to three 
dimensions which is clearly required for real ventricles. A naive extension of our mesh into 
a cubic array of sheets will, of course, succeed in achieving defibrillation. However, such 
an array of control sheets cannot be easily implanted in a ventricle. We have tried to see, 
therefore, if we can control the turbulence in a three-dimensional version of model (1) on 
a. L X L X Lz domain but with the control mesh present only on one L x L face. For the 
open faces we use open boundary conditions and for the other faces we use no-flux Neumann 
boundary conditions. Our defibrillation scheme works if < 4 (we have checked explicitly 
for L = 220) but not for > 4. A slight modification of this scheme is effective even 
for Lz > 4:-. Instead of applying a pulse for a duration Tc, we apply a sequence of n pulses 
separated by a time Tip and each of duration Tw We find that, if L = 256 and Lz — 8, 
defibrillation occurs in ~ 2002 ms with Tip = 22 ms, = 0.11 ms, n = 30 and a control 
pulse amplitude of 57.3 //A/cm^; if L = 128 and — 16, defibrillation occurs in ~ 1760 ms 
with Tip — 22 ms, Ty, — 0.11 ms, n — 15 and a control pulse amplitude of 57.3 //A/cm^. We 
further find that optimal defibrillation is obtained in our model if Tip is close to the absolute 
refractory period for model (1) without the Laplacian term. The efficacy of our control 
scheme in the three-dimensional case can be understood heuristically as follows: Control via 
a steady pulse docs not work for > 4 since the propagation of this pulse in the z direction 
(normal to our control mesh) is blocked once the medium in the interior of our simulation 
domain becomes refractory. However, if we use a sequence of short pulses separated by a 
time Tip, then, provided Tip is long enough for the medium to recover its excitability, the 
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control-pulse waves can propagate in the z direction and lead to successful defibrillation. 
Typical electrical defibrillation schemes use much higher voltages than in our study. Recent 
studies ll^ll^ have explored low-amplitude defibrillation methods in model systems. They 
constitute an advance over conventional methods, but lack some of the appealing features 



of our defibrillation scheme. For example, the scheme of Ref. |]T8[ works only when the 
slow variable (the analog of our g) is controlled; though this can be done, in principle, by 
pharmaceutical means, it is clearly less direct than control via electrical means. Reference 
||T9| uses electrical defibrillation, but has been demonstrated to prevent only one spiral from 
breaking up, as opposed to the suppression of a spatiotemporally chaotic state with broken 
spirals by our defibrillation scheme. We have checked explicitly that, for the spatiotemporally 
chaotic state of model (1), a straightforward implementation of the defibrillation scheme of 



Ref. [jT9| is ineffective. [This scheme applies pulses to the fast variable (e in model (1)) 
on a two-dimensional, discrete lattice of points.] The control current density in Ref. 
is comparable to our ~ 57.3/iA/cm^, which is much lower than 139/iA/cm^, the maximum 
value of the ionic current during depolarization in the Beeler-Reuter model. 
We have checked that our defibrillation scheme is not sensitively model dependent. For 
example, we have used the same scheme to eliminate spatiotemporal chaos associated with 
spiral breakup in the model for the catalysis of CO on Pt(llO) mentioned above We 
have found recently that our defibrillation scheme also works for the biologically realistic 
Beeler-Reuter model [0,|22| for VF. 

In conclusion, then, we have developed an efficient method for defibrillation by the elim- 
ination of spiral turbulence in model (1). Our method has the attractive feature that it 
uses very-low-amplitude pulses, applied only for a short duration on a coarse control mesh 
of lines. And, to the best of our knowledge, our study is the only one which shows how to 
attain defibrillation by the control of spatiotemporal chaos and spiral turbulence in a model 
for VF in both two and three dimensions. We hope our work will stimulate experimental 
tests of the efficacy of our defibrillation method. 
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Figure Legends 



Figure 1. Pseudograyscale plots showing the evolution of spatiotemporal chaos through 
spiral breakup in model (1) for physical times T between 880 ms and 1430 ms. The panels 
on the left show the excitability e and those on the right the recovery g for all points {x,y) on 
a two-dimensional L x L spatial grid with L = 128. The initial condition used is described 
in the text. 



Figure 2. Pseudograyscale plots of the e field in c? = 2 for L = 128 (top) and L = 
512 (bottom) illustrating defibrillation by our control of spiral breakup in model (1). The 
control mesh divides the domain into four equal squares for L = 128 and, for L = 512, 
into 49 equal squares, 28 equal rectangles along the edges and 4 small equal squares on the 
corners. For L = 128 we apply 57.3 yuA/cm^ from T = 891 ms to T = 935 ms and by T ~ 
1500 ms spatiotemporal chaos is all but eliminated (|e|, \g\ < 10^^^ at all grid points). For 
L = 512, we apply 250 /iA/cm^ from T = 55 ms to T = 759 ms to the e field. By T = 
1650 ms, spatiotemporal chaos is all but eliminated (|e|, \g\ < 10^^ at all grid points) (see 
http: / / theoryl .physics. iisc.ernet .in / rahul / images / heart .mpg] for an animated figure). 



Figure 3. Isosurface (e = 0.6) plots of initial states (left panels) and pseudograyscale plots 
of the e field on the top square face (right panels) illustrating defibrillation by the control 
of spiral breakup in model (1) for d = 3: (top) L = 128 and = 16, with the control mesh 
dividing the top face into 4 equal squares; and (bottom) L = 256 and Lz = 8, with the 
control mesh dividing the top face into 64 equal squares. In both cases we apply pulses on 
the control mesh [top face only) of 57.3 yuA/cm^ with Tip = 22 ms and = 0.11 ms (see 
text). In the former case, with n = 15, spatiotemporal chaos is all but eliminated by T = 
1760 ms (|e|, \g\ < 10~^ at all grid points); in the latter, with n = 30, it is all but eliminated 
by T = 2002 ms (|e|, l^f] < 10^^ at all grid points). 
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